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It is shown that temperature acts to disrupt the magnetization of Peierls distorted quasi-one- 
dimensional materials (QIDM). The mean- field finite temperature phase diagram for the field theory 
model employed is obtained by considering both homogeneous and inhomogeneous condensates. The 
tricritical points of the second order transition lines of the gap parameter and magnetization are 
explicitly calculated. It is also shown that in the absence of an external static magnetic field the 
magnetization is always zero, at any temperature. As expected, temperature does not induce any 
magnetization effect on Peierls distorted QIDM. 
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I. INTRODUCTION 

The recent experimental observation of a signature of one-dimensional (ID) transport in 30 nm wide graphene 
ribbons [l| certainly motivated the theoretical investigation of ID systems. A very interesting class of QIDM are 
the ones where a lattice distortion develops spontaneously with the appearance of a gap Aq, rendering the material 
an insulator. This is the mechanism of the Peierls transition^ which has been observed experimentally in various 
QIDM. Probably, one of the most well-known QIDM that exhibits this phenomenon is trans-polyacetylene (TPA). 
This polymer is a ID chain of CH groups with alternating single and double bonds, having one electron per site. In 
the tight binding approximation, TPA would be a metal. However, the interaction of the electrons with the lattice 
(also known as the spontaneous Peierls dimerization) is such that the energy gain in the system is always larger that 
the energy investment for distorting the lattice. As a consequence, polyacetylene is an insulator. After dopiiig, and 
the consequent vanishing of the gap at an experimentally measured critical doping concentration Uc = 6% [2| , the 
conductivity of TPA is enormously increased, presenting metal-like properties [3| . 

Although we shall refer to trans— GHx throughout the paper, due to the large amount of data on this material, 
the main results obtained here are robust and can be applied to any QIDM possessing the same dimerized structure 
of trans— CHx for which a field theory model is suitable to be employed [j, [^. The electron-phonon interaction 
in trans— CHx is described by the discrete Su-Schrieffer-Heeger (SSH) Hamiltonian and its continuum version, 
the Takayama-Lin-Liu-Maki (TLM) Hamiltonian Q. The TLM model is a relativistic field theory with two- flavor 
Dirac fermions that, with a doping equal to or greater than j/c, have linear dispersion relations with a Fermi velocity 
vp « lO^m/s, which is of the order of the velocity of the Dirac fermions in graphene Q. Employing the Gross- Neveu 
(GN) model |g|, that can be properly identified with the TLM model, as an effective field theory model for describing 
the insulator-metal phase transition in polyacetylene [lol - [l2j . there has been found a very good agreement with the 
experimental value, yc = a/ie = 6% [HI, where TV (=2 for TPA) is the number of spin degrees of freedom of 
the (delocalized) tt electrons, h is the Planck's constant divided by 2tt, vf = kpfi/m, kp is the Fermi wavenumber, 
m = h? /2toa?, and a (= 1.22A for TPA) is the lattice (equilibrium) spacing between the x coordinates of successive 
CH radicals in the undimerized structure, and ~ ^ (Aq ~ 0.7eV for TPA) is the critical chemical potential at 
which the GN model undergoes a first-order phase transition to a symmetry restored (zero gap) phase p^ . 

Regarding our choice of using a continuum model, a few comments arc pertinent here. The TLM model can be 
derived from the SSH model by expanding its Hamiltonian about the Fermi surface and keeping terms only to lowest 
order in a/^, where ^ = fivp / ^ is the electronic correlation length. This is fulfilled for the materials we are interested 
in and will be described by the TLM (GN) model, for which ^ >> a. Besides, the clear advantages of employing 
continuum models such as the TLM (or GN) are the analytical solutions that they provide and the fact that field 
theory methods are suitable for the calculation of the effective potential, which is appropriate for the analysis of the 
phase structure of the model, as in this work. On the other hand, the disadvantages are that continuum models with 
relativistic dispersion relations have an electronic spectrum that is unbounded below, and that the acoustic modes 
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are lost in the continuum limit. The firstproblem is resolved by adopting a certain energy cutoff, and the second, if 
only terms to next order in a/^ are kept [3|. 

In the zero temperature phase diagram of ID TPA under asymmetric doping, defined as an imbalance between 
the chemical potentials of the electrons with the two possible spin orientations ( "up" =t, and "down" =1) introduced 
in the system by the doping process, has been studied. As emphasized in |14| . the chemical potentials asymmetry 
between the f and ^ electrons can be achieved experimentally by the actuation of an external static magnetic field Bq 
on the polyacetylene wire, which breaks the spin-1/2 SU (2) symmetry. In the continuous model that describes 
the electron-phonon interactions in TPA has been introduced and the magnetization, the critical magnetic field Bq c 
at which there is a quantum phase transition to a fully polarized (magnetized) phase, and the magnetic susceptibility 
at zero temperature, within the field theory approach, have been obtained. 

In this paper, we study the thermal effects on the magnetic properties of Peierls distorted QIDM and verify the 
possibility of the existence of this fully polarized phase at finite temperature. The mean-field finite temperature phase 
diagram for the field theory model employed is obtained. The tricritical points of the second-order transition curves of 
the gap parameter and magnetization are explicitly calculated by considering both homogeneous and inhomogeneous 
A(a;) condensates. One of the main results of this paper is the demonstration of the "stationarity" of the tricritical 
point of the second-order transition line of the gap parameter under the the influence of an external (constant) Zeeman 
magnetic field. In other words, in a Peierls distorted QIDM under the influence of an external Zeeman magnetic field, 
the tricritical point obtained considering homogeneous condensates remains at the same location when inhomogeneous 
A(a;) condensates are taken into account. 

The paper is organized as follows. In Section II we introduce the model Lagrangian describing polyacetylene. In 
Section III the temperature dependent renormalized effective potential is presented. In this section we obtain an 
analytical expression for the effective potential at high temperature, as well as the chemical potentials dependent gap 
equation and the critical temperature at which the gap vanishes. Besides this, the tricritical points of the second- 
order transition lines of the gap parameter and magnetization are explicitly calculated. The temperature dependent 
magnetic properties of Peierls distorted QIDM are also obtained in this section. In the Summary we present the 
conclusions. 



II. MODEL LAGRANGIAN 

For the benefit of the reader, let us reproduce from the model Lagrangian and the basic definitions necessary 
to describe polyacetylene and equivalent Peierls distorted QIDM. As we mentioned already, the electron-phonon 
interaction in TPA is represented by the SSH model 0, which has a continuum version known as the TLM model 0. 
The TLM Lagrangian density in the adiabatic approximation (considering static configurations for which dA/dt = 0) 
is given by 

N , 

£tlm = V {ihdt - ihvFi^d, - 7oA(.T)) ^1 - -—- A2(x) , (1) 

^ 27rftz;FATLM 



where ?/; is a two component Dirac spinor = I ^ , representing the "left moving" and "right moving" electrons 

close to their Fermi energy, respectively, and j is an internal symmetry index (spin) that determines the effective 
degeneracy of the fermions. We define 1 =t, and 2 =4,. The gamma matrices are given in terms of the Pauli matrices, 
as 7o = cTi, 75 = —173, and A(a;) is a (real) gap related to lattice vibrations. Atlm = ^^j^ is a dimensionless 
coupling, where a is the 7r-electron-phonon coupling constant of the original SSH Hamiltonian, and K is the elastic 
chain deformation constant. The equivalence between the TLM and the Gross-Neveu (GN) model @, is established 
by setting Atlm ~ iv?- Note that the electron-phonon interaction term in Eq. ([T|) is an analog of the fermion- 
boson interaction in the field theory context, which appears in different models and dimensions. In four space-time 
dimensions, for example, this interaction has been investigated in the framework of the linear sigma model at finite 
temperature [T5| . 

The GN model has been investigated earlier at finite temperature and density several times (see for instance 

Refs. [ilHIiS), 

and recently considering also finite corrections to the leading order in the large N approximation [l^ . 
However, these approximations did not consider the effects of an external Zeeman magnetic field applied on the 
system, which is of fundamental importance in many physical situations, as in the investigation of metal-insulator 
transitions [l^ and magnetization in 2D electron systems [T9| . 

In order to consider the application of an external Zeeman magnetic field to the system and its effects, it is convenient 
to start by writing the grand canonical partition function associated with £gn, 
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Z= T>ij Tup exp { dr dx [Lqn 



(2) 



where ip = ip^'^o, P = l/fc^T, ks is the Bohzniann constant, and Xgn is the Euchdean GN Lagrangian density: 



Lgn = ^ ipj[-jo^dr + ifiVFlA - A(a;) + 7oAij]''/'j 



1 



A^x), 



(3) 



where jj,^ = fl + 5fi, /ij, = /l — (5/i. The Zeeman sphtting energy is given by AE = S^gy^sBQ [20| . where S'^ = ±1/2 
is the effective g-factor and /i^ = eTi/2ra 5.788 x 10^^ cV T^^ is the Bohr magneton, giving = ^iibBq. In [1 
we also have chosen /2 = /ic. 

Integrating over the fermion fields leads to 



Z = exp 4 / A^{x)\ U^.j^detDj, 

[ nVFAGN J J 



(4) 



where Dj = ^"fodr + iT^-vplidx + 7oMi ~ ^(2^) is the Dirac operator at finite temperature and density. Since A(a;) 
is static, we can transform Dj to the w„ plane, where aj„ = (2n + l)7rT are the Matsubara frequencies for fermions, 
yielding Dj = {—iuJn + Mj)7o + i^'^Fli^x ~ A(a;). After using an elementary identity \n{det{Dj)) = Trln(£)j), one 
can define the bare effective action for the static A(x) condensate 



/3 



2 

/ dx A^{x)+^i:v\n{Dj) 

•' 3 = 1 



(5) 



where the trace is to be taken over both Dirac and functional indices. The condition to find the stationary points of 
Seff[A] reads 



5S, 



effl 



SA{x) 



= = 



2/3 



GN 



-A{x) 



SAix) 



J^TvlniD,) 



(6) 



The equation above is a complicated and gene rally unknown functional equation for A(a;), whose solution has been 
investigated at various times in the literature j2ll426| . Its solution is not only of academic interest, but has direct 
application in condensed matter physics as, for example, in jsl, I25I], and in the present work. 



III. THE RENORMALIZED EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 

A. Homogeneous A(a;) Condensates 



For a constant A field the Dirac operator reads Dj — (— + Mj)7n + iTivp^iP — A, so the trace in Eq. ([5]) can 
be evaluated in a closed form for the asymmetrical 7^ 0) system [ij]. From Eq. (jH) one obtains the "effective" 
potential Vg/ f — ^ In 2., where L is the length of the system: 



K//(A,/.t,;,T) 



-A^ 



F^GN 



+ CXD 



dp 
2Trh 



2pEp + In ( 1 + er/^^r ) + In ( 1 + e^'^^t 



(7) 



Ind + e-'^^n+lnfl + e-'^^' 



where Ef .^^ = Ep± Ep = y/vj,p^ + A'^. 

The first term in the integration in p, corresponding to the vacuum part {fJ.^,i = T = 0), is divergent. Introducing 
a momentum cutoff A to regulate this part of Vef / , we obtain, after rcnormalization, a finite effective potential 
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K//(A) 



hve V A, 



3 

2^ 



ttTivf ymp 



(8) 



where is an arbitrary renormalization scale, with dimension of energy. The minimization of Ve//(A) with respect 
to A gives the well-known result for the non-trivial gap Q: 



Aq = mpe 



(9) 



From this gap equation we see that with the experimentally measured Aq and a, to and K which enters Atlm, one 
sets the value of mp. Equation ([SJ can be expressed in a more convenient form in terms of Aq as 



VeffiA) 



2TThvp 







Kl) 


- 1 





(10) 



which is clearly symmetric under A —> —A, which generates the discrete chiral symmetry of the GN model. As has 
been pointed out before 2l| , this discrete symmetry is dynamically broken by the non-pcrturbative vacuum, and thus 
there is a kink solution interpolating between the two degenerate minima A = ±Ao of ()10|) at x = ±oo: 



A(a;) = Aotanh(Aox). 

In the next subsection we discuss the effects of space dependent A(a;) condensates. 
We can rewrite Ve//(A, T) as 



(11) 



where 



dp 





In f 1 + e^'^^i' 1 + In f 1 + e^'^^t ) + ln(l + + In f 1 + e 



(12) 



(13) 



Since we can not calculate expression ()13p in a closed form, we shall use a high temperature expansion to evaluate it. 
Using the function 



I{a,b) 



dx 



In 1 + e 



In 1 + e 



(14) 



where a = A/ksT, and b = ^/kBT, which can be expanded in the high temperature limit, a << 1 and b << 1, 
yielding, up to order a* and b^ ^t\ . 



/(««M«l) = y + y-yln©-^(l-7.)- 



186^(5) 



1287r4 



fe^a^ + {a%^) , (15) 



where jp ~ 0.577... is the Euler constant and ^(n) is the Riemann zeta function, having the values ^(3) ~ 1.202, and 
^(5) « 1.037. With the equation above, together with Eq. ([SJ, the high temperature asymmetrical effective potential 
is written as 



1 



A2 



nhvf 



In 



TrksT 

"aT 



3hvj 



-{keTf 



(16) 



2'K'hvi 



7g(3) A^ _ 7^(3) 
87r2 {keTY " 47r2 



/ 2 , 2\ ^ 



186^(5) 
647r4 



A^ 



The equation above may be rearranged in the form of a Ginzburg-Landau (GL) expansion of the grand potential 
density, which is appropriate to the analysis of the phase diagram in the region near the tricritical point, 
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Veff = ao + a2A^ + a4^A'^ 



(17) 



where 



a2(m>^) 
a4(m'^) 



(18) 



1 



TlflVF 



In 



^fcsTX ^ 7^(3) (A^t + /^i) 



-14^(3) 



186^(5) (Mt + m') 

47r2 (fcBT)2 



(19) 



Extremizing Vg// we find the trivial solution (A = 0) and the chemical potential and temperature dependent gap 
equation 



a2 



(20) 



which has meaning only if the ratio ^ is negative. Besides, a stable configuration (i.e., bounded from below) requires, 
up to this order, > 0. At the minimum Veff reads 



ef f,min — ^^0 



4a4' 



(21) 



The critical temperature Tc is, by definition, the temperature at which the gap vanishes. Thus, at Tc we have 
a2 = or 



In 



ef-EAn 



7^(3) (/^? + f^l) 

8^2 {kBT,)^ 



0, 



(22) 



where Tc = Tc{fi^, /ij.). As will become clear below, the equation above defines a second-order transition line separating 
the non-metallic (A 7^ 0) and metallic phases (A = 0). At /j,^ = /j,^ = 0, we recover the well-known result for the 
temperature at which the discrete chiral symmetry is restored psl ]: 



e'''^ An 

TT kB 



(23) 



In order to find Tc(^-f-,^^) we define dimensionless variables 77 
Eq. ([23]) we rewrite the L.H.S. of Eq.(l22l) as 



7^(3) (K+MI) , , 

871-2 (/cj3Tc(0))2 "^^^^ 



^ , and with the help of 



7^ 



yit) = ln(i) 



(24) 



The zeros of y{t) for a given rj, i.e., for a given fi^ + jj,^, are the respective Tc. This defines the (second-order) Tc versus 
fi^ + /xj phase diagram. As can be seen in Fig. [U the graphical analysis of y{t) shows that there is no solution for this 
function for t] above certain value, that we define rjtc- Besides, at Tjtc we have y = y' = 0. These two equations give 
ttc and rjtc for the tricritical point Ptc = {iltcittc)- Strictly speaking, y and y' are associated with the coefficients of 
the second-order and forth-order terms of the effective potential expanded in powers of A [2^ . Solving the equations 
y ~ Q and t/' = sclf-consistcntly (which is equivalent to solve a2 — on — Q) ^ we obtain 



7g(3) (Mt + ^^l)tc 
8^2 ikBTciOW 



1 

2e 



tt 



1 



(25) 
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For r] above certain value and less than rjtc, the function y presents two solutions (not shown in Fig.[T]) for Tc- However, 
the lower of these always corresponds to unstable solutions. The second-order transition curve, defined as the line 
starting at the point (0,Tc(0)) and ending at the point + /i|)(c,Ttc), comes simply from the solution of the 
gap equation. This curve, shown in Fig. [5J represents a system at finite temperature where the chemical potentials 
/i| + /^l = 2(/2^ + S^^) start from zero and increases until + /ij)(c. Note that p,^ + S^j? is zero if and only if /I^ 

and 5y?' arc both zero. It is well known that below the tricritical point one has to properly minimize the effective 
potential rather than using the gap equation as the transition becomes first order. Thus Eq. (j22[) cannot be used for 
finding below Ptc since this equation is valid only for the second-order transition. In this case Tc has to be find 
numerically, through the equality Veff{n-f, /ij., A = Amin, Tc) = VeffifJ-^, fJ-i, A = 0, Tc), where A,„i„ is the non-trivial 
minimum of Vef / . 




FIG. 1: The function y{t), as a function of t = jrj^, for different values of rj. The bottom curve is for 77 = 0.01, the second 
curve is for rj — rite ~ ~ 0.184, and the top curve (with no solution) is for 77 — 0.4. 

The number densities n^^^ — — ^ Vef / (A, /^-f,^, T) read 

dp 



where nk{E^"^ ) 



"t,4. +1 
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is the Fermi distribution function. The density difference 



(26) 



5n = — 



(27) 



is zero if = ^/.igBo = 0, at any temperature, since in this case we have the equalities nk{E^) = nk{E^), and 
nk{E^) = nk{E^). The physical meaning of these results is that at zero external Zeeman magnetic field, the t (up) 
and J, (down) electrons of the conduction (+) band have the same density, and the same for the electrons of the 
valence (— ) band. 

In the high temperature regime, the number densities are given by 



nt:liT) 



ttTivf 



7^(3) A2 



47r2 (fcBT)2 



(28) 
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FIG. 2: The phase diagram t = tToJ ^ ^ function of 77 = "^^^^J (^kB^Tio))'-^ ' ^^^^ El"?- (ESJ- The smaU dot at the end of the 
second order transition line represents the tricritical point Ptc = ■^). Below this point the transition is of first order. 



In the high temperature limit the total number density, nT(T) = n^{T) + nj^(T), is independent of the applied field 



nhvp 



7^(3) A2 



47r2 {keTy 



(29) 



and for the density difference we obtain 



xflVf 



1 



7e(3) A2 



47r2 [keTY 



(30) 



that, as we have observed before, is clearly zero if Bq ~ Sfi = 0. Since the densities have to be evaluated at the 
minimum of the effective potential, we use Eq. pO| in the equation above and find the temperatures at which the 
densities and, consequently, the density difference vanish. These temperatures arc the solutions of 



IE 



In 



Aq \ 7g(3) (/i^ + M^) 
^kBT*) 8^2 {kBT*Y 



= 0, 



(31) 



where T* = T* {fi-f , fi^) . As for the gap parameter, the equation above defines the second-order line for the densities 
and the density imbalance. At fi-^ ~ fi_i = 0, we get 



7^c(A^t = Mi = o)^r;(o) 



-2 Ao 



(32) 



It is very easy to see that T*{0) 
before we find 



Te(0) 



which coincides with Ttc, where Ttc is given by Eq. (|25l). Proceeding as 



7^(3) + ^^l)tc 

87r2 (fcBr*(0))2 



1 

2~2 



^ o 1 tin 



T*(0) 



(33) 
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defining the tricritical point for the densities, total density and density imbalance second order curves. 

Let us now verify the possibility of a fully polarized state at finite temperature. It would be possible with a magnetic 
field with a intensity such that in Eq. ([28|) vanishes. In this case = fi — Sfic = /^c ^ f A'-B-Bq^c — 0, or 



Bo.c^^, (34) 



yielding, for TPA (for which fx^ ^ ^ and 5 « 2) a critical magnetic field 



Bo,c « 8.6 kT, (35) 

which is, as the critical magnetic field found in [l^ at zero temperature (~ 4.6 kT), a magnetic field of very high 
magnitude, compared to the maximum current laboratory values [sot . 



B. Inhomogeneous A(2;) Condensates 

Since we consider the addition of a chemical potential (i.e., doping) in the theory representing Peierls distorted 
QIDM, and the effects of a Zeeman magnetic field on these materials, some important remarks are in order. It 
is well-known that doping in conducting polymers with degenerate ground states results in lattice deformation, or 
non-linear excitations, such as kink solitons and polarons, meaning that A(a;) can vary in space [sl. [Sll. [3^ . Therefore, 
one may expect not only homogeneous- like configurations (as considered in the previous subsection), but also that 
the inclusion of these excitations in any theoretical calculation in this model should be considered. In this context, 
within the GN field theory model that we are considering, by taking into account kink-like configurations in the large 

approximation, the authors of Refs. [23l - |26j found evidence for a crystalline phase that shows up in the extreme 
T ~ and large /i part of the phase diagram, while the other extreme of the phase diagram, for large T and small ^, 
seemed to remain identical to the usual large N results for the critical temperature and tricritical points, which are 
well-known results [l^ for the GN model. 

To take into account the effects of inhomogeneous configurations in the GL expansion of the grand potential density, 
let us write Eq. ([T7]) in terms of A(x) and its derivatives up to [2^-[2^: 

Veffix) = ao + a2A(x)2 + a4[A{x)^ + Aix)'^], (36) 

where A(a:;)' = dA{x)/dx. A straightforward variational calculation gives the following condition for the minimization 
of the free energy E = J Veff{x) dx: 

A(a;)" - 2A(a;)3 - ^A(x) = 0. (37) 
The general solution of an equation of the form 



A(a;)" - 2A{xf + (1 + iy)AlA{x) = 0, (38) 

can be written as 



A(x) = Ao\^sn(Ao2:; v), (39) 

where sn is the Jacobi elliptic function with the real elliptic parameter < v < 1- The sn function has period 2K(i'), 

where K(z^) = /J^^^[l — iysin'^{t)]^^/'^dt is the complete elliptic integral of first kind. A{x) in ((39)) represents an array 
of real kinks. When ^ = 1 Eq. is reduced to the single kink condensate given in Eq. pT|) . By comparing Eqs. ([57)1 
and (|38p one can identify the scale parameter Aq as 



(40) 
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Given that > 0, the solution for inhomogeneous condensates has physical meaning only if the ratio ^ is negative, 
as in the case of homogeneous condensates. In terms of Eq. (|39|) it can be shown that the x dependent grand potential 
density can be written as 



Veffix) = ao + a2A{x)^ + 04- [(1 + iy)AlA{x)^ + i^A^] . 



(41) 



Averaging over one period, it is found [2^ that < A(a;)^ >= (l — km I ^O' where E(i/) is the complete elliptic 



integral of second kind. The ratio E(i/)/K(z/) is a smooth function of 1/ interpolating monotonically between and 1. 
Thus wc can write 



< Veffix) >= ao + A2AI + A4A 



where 



0' 



(42) 



A2 

A4 



a2 I 1 
1 



EM 
K(^) 

f (1 + 



(43) 



EH 
K(^) 



The interesting results obtained considering inhomogeneous condensates are: 



(1.) For — 1, ^l^^^j — 0, so the grand potential density is that of the homogeneous case, Eq. ((T7)) . at the 
non-trivial minimum. 

(2.) For I' = 0, = 1, so the grand potential density is that of the metallic phase, for which A = and 
Veff = ao- 

(3.) The tricritical point is still found for a2 = 0:4 = 0. These coefhcients are /i-j-.j, and T dependent and were not 
affected by the space dependence of the condensate A{x). Then the location of the tricritical point in a Peierls 
distorted QIDM under the influence of an external Zeeman magnetic field is unaltered even considering a x dependent 
grand potential density. This happens because in the high temperature limit the influence of the Zeeman field is not 
sufficient to change the position of the tricritical point. This same conclusion has been obtained for the symmetric 
case (Sfi = Bo^O) U^-^M- 



C. Magnetic Properties 

The Pauli magnetization of the chain in the high temperature limit has the following expression: 



Mhigh t{T) = ^iBSUhigh t{T) 



1 



7^(3) A2 



47r2 {ksT)'^ 



In 



T 



T*(0) 



7^(3) (a^? + 



Bn 



where we have made use of Eq. (j20p to leading order in -j^^^^. The second-order line where the magnetization 
vanishes is the same as the one given by Eq. (j3ip . Finally, we obtain the magnetic susceptibility in this regime 



(44) 



Xhigh t{T) — 



dM^gh t{T) 
dBo 



= X(0)+X(T), 



(45) 



where 



(46) 



10 



and 



X{T) 



nhvF 



In 



T 



7^(3) 



Tc{0)J 47r2(fcBT)2 



2 2 o2 



(47) 



x(0) is the well known zero temperature contribution for the Pauli expression of the magnetic susceptibility for nonin- 
teracting electrons. The function Xhigh t {T) also behaves at finite temperature as the densities and the magnetization, 
with a second-order transition up to a tricritical point given by Eq. (j33[) . Below this point the transition is again of 
first order. 

Note that in spite of the fully polarization, at Bq c the magnetization is given by exactly the same expression shown 
in Eq. dM]). With the help of Eq. ^ we find: 



Mhigh tAT) = l-J-sn^iT) 



In 



7^(3) A2 1 ^ 
T \ 7^(3) (m^ + 



T*{0) 



(48) 



This shows that this function is indeed continuous for < Bq < Bq ^ 



IV. SUMMARY 



We have investigated the mean-field finite temperature phase diagram of QIDM under the influence of an external 
Zeeman magnetic field. We found that the gap parameter and the magnetization (as well as the densities and 
density imbalance, and the magnetic susceptibility) of asymmetrically doped QIDM have a similar second-order 
behavior until their respective tricritical points are reached. Below these points the transitions are of first order. 
We found these two tricritical points analytically. He have shown that the location of the tricritical point in the 

t j^jj^^ versus 77 ''1^^'' (fc^T^^o))^ ) phase diagram stays at the same place by considering both homogeneous 
and inhomogcneous condensates, as occur in symmetric = Bq = systems [23l - [26j . We have shown that for 
the particular case of TPA, in order to have a fully polarized organic conductor at finite temperature it would be 
necessary to have a very high critical magnetic field, namely Bq c- However, for a given magnetic field below i?o.cj 
partial polarizations (magnetizations) can be realized experimentally, provided the temperatures are kept outside the 
"non — metallic^' region of Fig. [21 It is worth noting that, according to Eq. for other ID systems with a smaller 
critical chemical potential or with a greater effective g- factor, a smaller (attainable) critical magnetic field necessary 
for a fully polarization of the QIDM would be found. As a final remark, it would also be very interesting to study the 
transport properties of the asymmetrically doped Peierls distorted QIDM at zero and finite temperature, employing 
the field theory approach. We intent to address these topics elsewhere. 
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